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Abstract 

Using ab-initio calculations based on density functional theory, we investigate the effects of 
vacancies on the electronic and magnetic properties of zigzag SiC nanoribbons (Z-SiCNR). Single 
(Vc and Vsi) and double (VsiVsi and VgiVc) vacancies are observed to induce magnetism in Z- 
SiCNRs. The presence of a single Vsi does not affect the half-metallic behavior of pristine Z-SiCNRs; 
however, a single Vc leads to a transition from half-metallic to metallic behavior in Z-SiCNRs due 
to the edge Si p orbitals and the atoms surrounding the vacancy. The interactions of vacancies with 
foreign impurity atoms (B and N) are also investigated and it is observed that VsjiVc does not only 
suppress the oscillatory type magnetism of VsjVc, but also retains the half- metallic character of the 
pristine Z-SiCNRs. The defect formation energies of vacancies can be reduced by substitutional B 
and N atoms. We believe that ferromagnetism is expected if Z-SiCNR are grown under suitable 
conditions. 

PACS numbers: 73.22.-f, 71.55.-i, 75.75.-c, 71.15.Mb, 71.15.Nc 
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I. INTRODUCTION 



Silicon carbide (SiC) is an attractive material for numerous technological applications, 
mainly in harsh environments. Bulk SiC is known to possess outstanding properties (such as 
high thermal conductivity, high breakdown electric field, high electronic mobility, excellent 
chemical and physical stability, good radiation resistance, and wide band gap) 1,2 which make 
it a suitable semiconductor for high-power, high-temperature, and high-frequency devices. 
Furthermore, SiC nanowires and nanotubes, which have already been synthesized,- 1 ^ exhibit 
excellent characteristics and are good candidates for applications ranging from hydrogen 
storage medial and gas sensors^ to optical^ and field-emission devices.- 

In the last few years, the successful synthesis of SiC nanotubes^ and, recently, the theoreti- 
cal prediction of the stability of two-dimensional SiC monolayer with honeycomb structure^ 1 ^ 
have stimulated increasing interest in SiC nanosheets and nanoribbons (NRs). According to 
recent theoretical studies, SiC nanosheets as well as armchair SiCNRs (A-SiCNRs) behave 
as nonmagnetic wide band gap semiconductors,- - — whereas zigzag SiCNRs (Z-SiCNRs) are 
magnetic (with very small magnetic moments)^ 1 ^ and can present metallic or semiconduct- 
ing character, depending on the width of the nanoribbonii"— (half-metallicity was predicted 
for Z-SiCNRs narrower than 4 nm,- which makes these NRs promising candidates for spin- 
tronic applications). However, ab initio investigations have recently shown that these char- 
acteristics can be modified by the presence of some impurities and defects. For example, 
it has been observed that (i) substitutional B, N, As, and P impurities induce magnetism 
in SiC sheets,— (ii) half-metallic Z-SiCNRs become metallic when doped with N atoms,— 
and (iii) B (N) substituting an edge Si (C) atom transforms semiconducting Z-SiCNRs into 
half-metallic systems.— In addition, Si vacancy has been shown to induce magnetism in 
nonmagnetic SiC sheets^ 1 ^ and A-SiCNRs.— 

Intrinsic defects, especially vacancies, have been suggested to be related to the origin 
of magnetism in SiC structures: besides the aforementioned studies of vacancies in SiC 



sheets and A-SiCNRs (Refs. |10| andll6l). there are experimental evidences that defects dom- 
inated by Si+C divacancies can induce room-temperature ferromagnetism in diamagnetic 
SiC crystals.— 1 ^ Despite the great potential of Z-SiCNRs for spintronics, the role of intrinsic 
defects (vacancies) and their interactions with foreign impurities in these NRs have not, to 
our knowledge, been investigated so far. In this work we investigated, by means of ab initio 



calculations, the effects of vacancies on the electronic and magnetic properties of Z-SiCNRs, 
and the interactions of vacancies with foreign impurity atoms. Our extensive results indicate 
that the presence of one C (Vc) or Si (Vsi) vacancy per supercell, as well as two vacancies 
(VsiVsi or VgiVc), can induce large magnetic moments in Z-SiCNRs. The defect formation 
energy of Vs; is decreased when the interactions with B and N impurities are considered. 

II. COMPUTATIONAL DETAILS 

Spin-polarized calculations were performed within the framework of the density func- 
tional theory— as implemented in the SIESTA code,-22 employing the generalized gradient 
approximation of Perdew, Burke, and Ernzerhof^i for the exchange-correlation functional 
and norm-conserving Troullier- Mart ins pseudopotentials^ to describe the electron-ion in- 
teractions. We used an energy cutoff of 200 Ry for the real-space mesh and a double-zeta 
basis set with polarization functions for all atoms. We considered Z-SiCNR with width 
W = 6 (6Z-SiCNR), as depicted in Fig. [TJ The 6Z-SiCNR was simulated within the su- 
percell approach with 84 atoms in the unit cell (including H atoms to passivate the edge 
dangling bonds) and vacuum regions of 17 A and 20 A along the y and z directions, respec- 
tively. All atomic positions were fully relaxed until the forces on each atom were smaller 
than 0.05 eV/A. The Brillouin zone was sampled using Monkhorst-Pack fc-point meshes of 
llxlxl and 121 x 1 x 1 for the total-energy and electronic-structure calculations, respec- 
tively. Using the Mulliken population analysis, the local magnetic moment was calculated 
by m = p^ — p^, where p^ (pr) represents the spin-up (-down) valence electrons of an atom. 

In this work, we investigated the presence of a single (one Vc or Vsi per supercell) and 
double (VsiVsi and VsiVc) vacancies in 6Z-SiCNR. The energetic stability of these systems 
was examined by calculating their defect formation energies Ef. We used 

E f = £[V x /6Z-SiCNR] - £[6Z-SiCNR] + /i X , (1) 

for single vacancy (X = Si or C), and 

E f = - {E[V Si V x /QZ-SiCNR} - £[6Z-SiCNR] + /i Si + p x } (2) 

for double V Si V x (n = 2 if X = Si, and n = 1 if X = C). Here, £[6Z-SiCNR] is the total 
energy of pristine 6Z-SiCNR, and £[Vx/6Z-SiCNR] and £[V/ S i"\4/6Z-SiCNR] are the total 
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FIG. 1. (Color online) Structural model of 6Z-SiCNR. Labels 1,2, 16 indicate the positions of 
the vacancies or impurities. 

energies of 6Z-SiCNR with single Vx and double VsiVx, respectively. The chemical potential 
of the atomic specie X is /ix- In equilibrium conditions, 

Msi + Mo — ^ u c k = + Mc ulk " Ai? / (SiC), 

where Ai?/(SiC) is the formation heat of SiC. In the Si-rich limit, /isi = fJ^ lk and He = 
MsiC 4 ~ /"s? 1 ^ whereas in C-rich conditions /ic = McT 11 * an d Msi = l^sfc ~~ Mc 1 ' 15 - We considered 
the 3C-SiC structure in the calculation of i^q. To check the quality of the pseudopotentials 
and computational parameters which are used in the present calculations, we calculated 
AH f (SiC) = 0.679 eV, which is in good agreement with the experimental value of 0.72 eV.— 

III. RESULTS AND DISCUSSION 

First, we examined the electronic and magnetic properties of pristine 6Z-SiCNR. The 
electronic band structure [Fig. EJ^a)] and the density of states [Fig. |3]^a)] show that 6Z- 
SiCNR exhibits half-metallic behavior: the spin-up channel is semiconducting (with a direct 
band gap of about 0.41 eV at the T point) whereas the spin-down channel is metallic. 
As can be seen in Fig. [3^a), the partially occupied spin-down electronic state is mainly 
composed of p orbitals of the edge C and Si atoms. The 6Z-SiCNR also has a very small 
magnetic moment, about 0.03/is per supercell. Our results are in agreement with a previous 
theoretical work 11 which reports magnetic moment of 0.023/xb per cell in 6Z-SiCNR, and 
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predicts half-metallicity in Z-SiCNRs narrower than 4 nm. The spin-density distribution 
of pristine 6Z-SiCNR [Fig. IHa)] shows that the magnetic moments are mainly localized on 
the edge atoms, and their orientations are parallel at each edge and antiparallel between 
the two edge atoms. Each edge C (Si) atom has a local magnetic moment of about 0.16/is 
(— 0.17 he), while the total magnetic moment of the other C (Si) atoms is about 0.29/1^ 
(— 0.20/is). The strong contribution of the edge C and Si atoms to the magnetization of 
6Z-SiCNR can also be seen in Fig. [3ta): the spin-up p states of the edge C atoms are 
completely occupied whereas the spin-down states are partially occupied; for the edge Si 
atoms, we note that their spin-down p states have a small occupation whereas their spin-up 
states are empty. 

Next, the effects of vacancies on the electronic and magnetic properties of 6Z-SiCNR were 
investigated. We considered initially one Vsi or Vc per supercell, which corresponds to a 
defect concentration of 2.78%. Vsi and Vq were created at the center of the nanoribbon, by 
removing Si and C atoms labeled 1 and 7 [see Fig. [T] , respectively. The presence of a single 
Vsi and Vc induces magnetic moments of 4.19 and 1.19/is per supercell, respectively. We 



(a) (b) (c) 




FIG. 2. (Color online) Electronic band structure of (a) pristine 6Z-SiCNR, (b) Vsi/6Z-SiCNR, and 
(c) Vc/6Z-SiCNR. The band structure of each system is presented in two ranges: \E— Ep\ < 2.0 eV 
(top panel) and \E — Ep\ < 0.05 eV (bottom panel). Solid red and dashed blue lines indicate spin- 
up and spin-down bands, respectively. The Fermi level is set to zero and indicated by the solid 
black line. 
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FIG. 3. (Color online) Density of states (DOS) of (a) pristine 6Z-SiCNR, (b) y Si /6Z-SiCNR, and 
( c ) Vc/6Z-SiCNR. In the middle [bottom] panels the shaded regions indicate the DOS of the C 
[Si] atoms; solid red lines show p orbitals of the edge C [Si] atoms; and dashed blue lines show p 
orbitals of the three C [Si] atoms around Si 1 [C ] in (a) and Vg- [Vq] in (b) [(c)] . The Fermi level is 
set to zero and indicated by the solid black line. We considered a Gaussian broadening of 0.1 eV 
for the DOS diagrams. 
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FIG. 4. (Color online) Contour plot of spin-density distribution of (a) pristine 6Z-SiCNR, (b) 
Vsi/6Z-SiCNR, and (c) Vc/6Z-SiCNR. For comparison purposes the spin-density isosurface of 
V c /6Z-SiCNR system is also shown in (d) at 8.0 x 10" 4 e/Bohr 3 . 



found that one Vsi per supercell leads to an increase, from 0.03 to 4.19//B, in the magnetic 
moment of 6Z-SiCNR, which is similar to that reported for SiC sheets and A-SiCNRs in 
previous ab initio studies where a single Vsi was found to induce a large magnetic moment 
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of 4.0^£ in nonmagnetic SiC sheets^ 1 ^ and A-SiCNRs.— On the other hand, Vc appears to 
play different roles in the magnetism of Z-SiCNRs, A-SiCNRs and SiC sheets: while we found 
a magnetic moment of 1.19/ig induced by a single Vc in 6Z-SiCNRs, recent theoretical works 
have reported that the presence of a single Vc does not give rise to any magnetic moment 
in SiC sheets^ 1 ^ and A-SiCNRs.— This suggests that 6Z-SiCNR, where magnetism can be 
tuned either by Si or C vacancies, can be a promising candidate for spintronics. 

FigureBJb) presents the distribution of spin density in 6Z-SiCNR with single Vsi (Vsi/6Z- 
SiCNR). We note that the spin density is mostly localized around the vacancy and on the 
edge C atoms: the local magnetic moment of the each edge C atom varies from 0.233 to 
0.287/is, while the three C atoms surrounding the Vsi have magnetic moments of 0.435, 0.865 
and 0.865/is- This can also be seen in Fig. Etb), which shows that the spin-up p states of 
the edge C atoms and of the C atoms around the Vsi are completely occupied whereas their 
spin-down states are partially occupied. The presence of a single Vsi destroys the complete 
occupation of p orbitals of the C atoms surrounding the vacancy, giving rise to empty spin- 
down states above the Fermi level. The local magnetic moment of each edge Si atom is about 
— 0.04/iB, which suggests that the C atoms and the Si atoms couple antiferromagnetically. 
Such types of magnetic coupling are also observed in the other defect driven magnetic 
systems.— ~— The smaller magnetization of the edge Si atoms in V"si/6Z-SiCNR, compared 
to that observed in pristine 6Z-SiCNR (where we found magnetic moments of about — 0.17//B 
for each edge Si atom), can also be seen in Fig. [3j by comparing Figs. Eta) and (b) we note 
that the occupation of spin-down p states of the edge Si atoms is larger in pristine 6Z-SiCNR 
[Fig. Eta)] than in Vsi/6Z-SiCNR [Fig. Etb)]- Such changes in the occupation of spin-down 
p states also decrease the exchange splitting of the edge Si atoms which result in small 
magnetic moments. On the other hand, the exchange splitting is large in the edge C atoms. 

In the Vc/6Z-SiCNR system we found that the magnetic moment comes mainly from the 
edge C atoms [see Figs. Etc), Htc), andHtd)]. Each edge C atom has magnetic moment of 
about 0.20yU#, whereas the Si atoms surrounding the Vc have magnetic moments of -0.04, 
-0.07 and — 0.07/ie. It is interesting to see that both the Si and C atoms are isoelectronic 
where one could expect the same magnetism in both Si and C vacancies. However, we 
observed a quite different nature of magnetism, e.g., atoms surrounding Vsi (Vc) have large 
(small) local positive (negative) magnetic moments. This behavior can be understood due 
to different nature of wave functions of the p orbitals. Usually, cation/anion vacancies 
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induce spin positive polarization on the surrounding atoms, but Vc/6Z-SiCNR shows the 
spin density of the Si atoms is larger in the spin down states. 

The electronic band structure [Fig. 12(b)] and the density of states [Fig. Mb)} of Vsi/6Z- 
SiCNR show that the presence of a single Vsi does not affect the half-metallic character of 
pristine 6Z-SiCNR; in this system, the spin-down channel is metallic, whereas the spin-up 
channel is semiconducting, with indirect band gap of 0.41 eV. In contrast, we found that a 
single Vq transforms the half-metallic 6Z-SiCNR into a metallic system [see Figs. Etc) and 
E(c)]. Note that the half-metallicity is mainly destroyed by the p orbitals of the edge atoms 
and the atoms surrounding the vacancy. Figure [3]^c) shows that the partially occupied spin- 
up electronic state consists mainly of p orbitals of the edge Si atoms and of the Si atoms 
surrounding the vacancy, with no contribution from p orbitals of the edge C atoms. This 
shows that the electronic structure can easily be tuned either by creating holes (simply 
based on electron counting), i.e., Vc, or by electron doping, i.e., N-doping.— We attribute 
this electronic phase transition to hole doping. 2S The electronic phase transition is mainly 
due to Si-p orbitals which are delocalized as compared with the C-p orbitals. The total 
densities of states in Fig. [Hlja-c) show that the population of density of states (within an 
energy range of ±0.5 eV) in spin-up states of electrons increases as we move from pristine to 
Vc system. This increment is mainly caused by the edge Si atoms and Si atoms surrounding 
the vacancy. 

The energetic stability of single Vsi and Vc in 6Z-SiCNR was examined by comparing 
Ef of the Vsi/6Z-SiCNR and Vc/6Z-SiCNR systems. Our results for Ef, obtained using 
Eq. (pQ), are listed in Table U As can be seen, the occurrence of a single Vc in Z-SiCNRs is 
energetically favored over single Vsi- This behavior is similar to that observed in bulk SiC, 29 
SiC sheets and SiC nanotubes,— where Vq has also been reported to be more stable than 
V Si . 

In order to study the interactions of Vsi and Vc, we also investigated the presence of two 
vacancies per supercell. Since Vsi appears to play a more important role on the magnetism 
of Z-SiCNR than Vc we considered configurations with (i) two Si vacancies, V^V^, and (ii) 
Si + C vacancies, Vg-Vc- The index z in V£ indicates the position of the X vacancy [see 
Fig. d]. 

The formation energies and the total magnetic moments per supercell for VgjVgj and 
VgjVc in 6Z-SiCNR are presented in Table fl] (the defect formation energies were calculated 
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using Eq. (J2])). These results show that the presence of double Si vacancies can give rise to 
large magnetic moments in Z-SiCNRs (for example, we found magnetic moments of about 
3.9/is in Vs-Vgj/6Z-SiCNR and Vg*Vg[/6Z-SiCNR). However, the formation of these magnetic 
VgYMgj/GZ-SiCNR systems is quite unlikely, since they have very high formation energies [see 
Table [I]. It is also noticeable the Vsi — Vsi interactions reduce the magnetic moments as 
compared with the single Vgi system. The formation of Vg^V^, which also induce magnetic 
moments in 6Z-SiCNRs, is energetically favorable when compared with a single Vsi and 
double VgjVgi. For instance, V^Vq, V^Vq 2 , and Vg-V^ 3 have formation energies smaller than 
2.8 eV, and induce magnetic moments of 1.34, -0.93 and -0.57/^b in 6Z-SiCNR, respectively 
[see Table [I] . It is worth mentioning here that recent experimental works have shown that 
defects created in 6H-SiC bulk by neutron or ion irradiations are mainly formed by V^Vc 
divacancies and induce room-temperature magnetism in diamagnetic SiC crystals.—^ Such 
experimental results illustrate that VsiVc are the stable intrinsic defects in 6H-SiC, and we 
also found that VsiVc types defects are easy to be formed in NR as compared with Vsi- 
Table H clearly shows that VsiVc divacancies have oscillatory type magnetism where the 
magnetic moments oscillate and strongly depend on the location of the C vacancy. This rise 
and fall of magnetic moments in VgiVc divacancies system is similar to bulk SiC crystal.— 
Comparing the defects, we note that this rise and fall of magnetic moments is suppressed 
by filling the C vacant site with N, which will be discussed in the following paragraphs. 

Figure [5] and [6] present the electronic band structures [Fig. [5] and the spin-density distri- 
butions [Fig. [6] for some representative magnetic systems with two vacancies per supercell, 
viz., V^ 4 j/6Z-SiCNR and V^Vg/GZ-SiCNR. As can be seen in Fi g-El both V^V^/6Z-SiCNR 
and Vg-T^f /6Z-SiCNR systems exhibit semiconducting character, with band gaps of 0.91 eV 
(0.15 eV) and 0.31 eV (0.04 eV) for the spin-up (-down) channels of V^Vgl/eZ-SiCNR and 
Vg L i V3/6Z-SiCNR, respectively. These results indicate that the presence of double vacancies 
(VsiVsi or V^Vc) can transform half-metallic Z-SiCNRs into semiconducting systems. As we 
increased the vacancy concentrations to ~ 5.56% (in double vacancy systems), half-metal to 
semiconductor phase transition is observed. This phase transition is due to holes which are 
localized on the dangling bonds around the vacancies but also have some contribution from 
the edge atoms. Usually, when a vacancy is created some states around the Fermi level, 
which are occupied in the pristine system, are removed and such removal of states creates 
holes around the Fermi energy which changes the electronic structure of the material.—^— 
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TABLE I. The calculated formation energies, in Si-rich and C-rich conditions, and total magnetic 
moments per supercell (M) for single (Vsi and Vc) and double (Vg-Vg; and Vg-Vp) vacancies in 
6Z-SiCNR. The index i in V£ indicates the position of the vacancy [see Fig. [T]. Single Vsi (Vc) is 
localized on site 1 (7). 

Formation energies (eV) 



M (n B ) Si-rich C-rich 



Single vacancy 

V S i 4.19 8.053 7.374 

V c 1.19 2.526 3.205 



two Si vacancies 

viiVi 



-0.02 
3.86 
3.88 
2.99 
3.81 



2.806 
7.836 
7.439 
7.999 
7.481 



2.127 
7.157 
6.760 
7.320 
6.802 



Si + C vacancies 



V&Pg 0.06 4.871 4.871 

V^Vg 1.34 2.789 2.789 

V^Vg 0.14 6.256 6.256 

V^ 1.91 8.321 8.321 

V^Vg 1 1.99 9.305 9.305 

V&Vg 2 -0.93 1.949 1.949 

-0.57 2.593 2.593 

P&Vg 4 4.03 10.112 10.112 

VsiVg 5 2.30 7.388 7.388 

F S |Fc 6 3.47 9.305 9.305 



In Fig. [2(a) we see some states just below the Fermi level; however, when a Si vacancy is 
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FIG. 5. (Color online) Electronic band structure of the (a) Vg-V^/6Z-SiCNR and (b) 
VgjT/Q /6Z-S1CNR systems. The band structure of each system is presented in two ranges: 
\E — Ep\ < 2.0 eV (top panel) and \E — Ep\ < 0.05 eV (bottom panel). Solid red and dashed 
blue lines indicate spin-up and spin-down bands, respectively. The Fermi level is set to zero and 
indicated by the solid black line. 

created, some missing states (near to X-point) can be seen [see Fig. EJ^b)]. Similar holes 
were also observed in SiC monolayer with a Si vacancy.— In 6Z-SiCNR, the electronic and 
magnetic structures are not only modified by the atoms surrounding a vacancy but also by 
the edge atoms^ (also see Fig. 4(d)). This makes 6Z-SiCNR different from other materials 
where magnetism is caused by vacancies. Our results teach us that the electronic band struc- 
ture of Z-SiCNRs can easily be engineered by Vsi or Vq. The distribution of spin density in 
V^V^/GZ-SiCNR [Fig. H^a)] shows that the magnetic moments are mostly localized on the 
C atoms surrounding the Si vacancies. The C atoms around have magnetic moments of 
0.23, 0.88 and 0.88/ie, while the C atoms around have magnetic moments of 1.06, 0.91 
and 0.91/xb- On the other hand, in V^Vc/6Z-SiCNR [Fig. Mb)} we observe that the C atoms 
surrounding Vgj move to form new C-C bonds; thus, the dangling bonds of these atoms are 
recombined leading to almost zero local magnetic moments. In this system the magnetic 
moments are mostly localized on edge C atoms: we found magnetic moments between 0.16 
and 0.21/zb for each edge C atom. 

As our studied systems induce magnetism due to vacancies, and the magnetic moments 
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FIG. 6. (Color online) Spin-density distribution in (a) V^V^/eZ-SiCNR and (b) V^Vg/6Z-SiCNR. 

are not carried by conventional magnetic elements (Fe, Ni, or Co) where one can also calcu- 
late the true ground magnetic state of a material. Either Vs; or Vq in 6Z-SiCNR distorts the 
bond lengths which can give different atomic magnetic moments at different atomic sites. 
To get a qualitative analysis of magnetism in 6Z-SiCNR, we also carried out calculations by 
considering the antiferromagnetic (AFM) and ferromagnetic (FM) types of Vsi — Vsi inter- 
actions. We found that the FM state is more stable than the AFM state by ~ 0.14 eV per 
cell which is much larger than room temperature. Therefore, room-temperature magnetism 
is expected in 6Z-SiCNR. The FM stability against AFM was also confirmed through fixed 
moment calculations. It is emphasized that such kind of magnetic coupling strongly depends 
on the distance between the vacancies. 

Among the defects investigated, single Vsi was found to induce the largest magnetic 
moment [see Table H], which suggests that this defect can play an important role on the 
magnetism of Z-SiCNR. However, the formation of single Vsi in Z-SiCNR is expected to be 
energetically unfavorable due to its high defect formation energy. To explore this system, 
we need to reduce its defect formation energy and for this reason the role of substitutional 
B and N impurities in 6Z-SiCNR is also investigated. Since Costa and Morbeo^ recently 
reported that B prefers to occupy a Si site whereas N preferentially substitutes a C atom, 
we examined configurations with single Vsi a t position 1 and (i) B at Si site i — 2,..., 6 [see 
Fig. Op, Vg^i/6Z-SiCNR, or (ii) N at C site % = 7, 16, M£iVj/6Z-SiCNR. The energetic 
stability of these systems was determined from their formation energies 

E f = E[^Y^6Z-SiCNR] - £[6Z-SiCNR] + fi Si + fi x - fiy, (3) 

where X=Si and Y=B for the V^B^/eZ-SiCNR systems, and X=C and Y=N for V^/ 6Z-SiCNR. 
We used a-Boron bulk and N 2 molecule to obtain B and N chemical potentials, respectively. 
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TABLE II. Calculated formation energies, in Si-rich and C-rich conditions, and total magnetic 
moments per supercell (M) for 6Z-SiCNR with a single Vsi and B (N) at Si (C) site. The index i 
in V^B^ and V^N^ indicates the substitutional Si or C site for B or N impurity [see Fig. [TJ. 

Formation energies (eV) 



M (hb) Si-rich C-rich 



B at Si site 

V^Bl 2.86 6.904 5.546 

V^Bg 4.90 8.235 6.877 

2.24 7.058 5.700 

V^B|i 4.61 8.611 7.253 

V^Bg 4.81 7.565 6.207 



N at C site 

V^Nl 0.63 3.240 4.871 

V^iVg 3.06 6.353 6.353 

V^JVg 3.07 5.769 6.256 

V&Ng 3.01 6.690 6.690 

V^Ng 2.38 5.295 5.295 

V&Ng 3.12 5.524 5.524 

V&Ng 3.12 5.789 5.789 

V^Ntf 3.09 5.767 5.767 

V&Ng 3.04 6.319 6.319 

V&Ng 2.95 7.063 7.063 



The formation energies (in Si-rich and C-rich conditions) and the total magnetic moments 
per supercell for V^^i/6Z-SiCNR and V^iV£/6Z-SiCNR systems are listed in Table M 
From these results, we note that some Vgj-Bgj/GZ-SiCNR systems, and all Vg^A^/GZ-SiCNR 
configurations are energetically more favorable than V^/GZ-SiCNR, which indicates that one 
B atom occupying a Si site, as well as one N atom occupying a C site, can stabilize a single 
V Si in Z-SiCNR. In addition, we observe that the V^./GZ-SiCNR and V^iV^/6Z-SiCNR 
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systems have large magnetic moments (except for VgjiVc/6Z-SiCNR whose magnetic moment 
is 0.63/^b). This smallest magnetic moment can be expected in Vg\iVc/6Z-SiCNR based 
on the separation between Vg- — Nq which is 1.81A. We found magnetic moments larger 
than 2{iB for all other configurations with Vgi-Bg; or VgiiVc- However, the VsiiVc/6Z-SiCNR 
systems are more interesting than Vs;-Bsi/6Z-SiCNR. When N is doped at the vacant C 
site in VsiVc/6Z-SiCNR, then the C-Vsi-types interactions increase the magnetization. So, 
we believe that N in 6Z-SiCNR has twofold roles; it not only stabilizes the intrinsic defects 
but also increases the magnetization of 6Z-SiCNR. From the whole thermodynamic data, 
it is summarized that those systems with low magnetic moments have smaller formation 
energies. 

The electronic band structures for the Vg^-Bgj/GZ-SiCNR and Vg^A^/GZ-SiCNR systems 
are depicted in Fig. As can be seen, the semiconducting character of V^V^/GZ-SiCNR 
is maintained when B is doped at the vacant Si 4 site: both the spin-up and spin-down 
channels of Vg^Sgj/GZ-SiCNR [Fig. [7(a)] are semiconducting with band gaps of 0.19 and 
0.13 eV, respectively. However, the presence of a N atom occupying a vacant C site can lead 
semiconducting VsiVc/6Z-SiCNR systems to turn into half- metallic ones. Fig. [7(b) shows 
that the V^iVc/GZ-SiCNR system is a half- metal: while its spin-down channel is metallic, 
the spin-up channel is semiconducting with band gap of 0.52 eV. It is interesting to note 
that Vg'i-Bgj/GZ-SiCNR remains semiconducting consisting with our picture of hole-induced 
electronic phase transition. When a vacant Vgj site is filled with boron in bivacancy Vg-Vgj 
system, it means that we are filling the empty orbitals (holes) created by Vgj. However, 
the Vj^-Bgj system has sufficient holes, compared with the pristine 6Z-SiCNRs, and remains 
semiconducting as expected. Boron mainly induces impurity states near the Fermi level in 
both spin-up and spin-down channels [see Fig. [7(a)] . On the other hand when N is doped at 
Vc site in V^Vq system, electrons are injected into V^Vq and we have sufficient electrons to 
fill the empty states. These electrons are mostly in the spin-down states and have sufficient 
energy to cross the Fermi energy, and V^Nq, e.g., transforms to half-metal as compared 
with VgjVQ. Such re-entrant behavior of 6Z-SiCNR is ascribed to electron doping. 

The spin-density maps presented in Fig. [8] show that the presence of a B (N) impurity fill- 
ing a vacant Si (C) site strongly affects the distribution of spin-density in VsiVs;/6Z-SiCNR 
(VsiVc/6Z-SiCNR). When a B atom occupies the vacant Si 4 site [see Fig. [8(a)] the dan- 
gling bonds of the C atoms surrounding Vgj are saturated leading the magnetic moments of 
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(a) (b) 




FIG. 7. (Color online) Electronic band structure of the (a) Vg|.Bgj/6Z-SiCNR and (b) 
VgiiVc/6Z-SiCNR systems. The band structure of each system is presented in two ranges: 
\E — Ep\ < 2.0 eV (top panel) and \E — Ep\ < 0.1 eV (bottom panel). Solid red and dashed 
blue lines indicate spin-up and spin-down bands, respectively. The Fermi level is set to zero and 
indicated by the solid black line. 




FIG. 8. (Color online) Spin-density distribution in (a) V^I/GZ-SiCNR and (b) V^iVg/GZ-SiCNR. 

these C atoms to decrease from 0.23, 0.88 and 0.88/i B in V£V$/6Z-SiCNR to 0.004, 0.039 
and 0.039//B in yg^gj/GZ-SiCNR. In addition, two of the three C atoms surrounding Vg\ 
move closer to each other to form a C-C bond, which leads to a decrease (from 0.91/ie in 
Vg^Vg^/6Z-SiCNR to 0.02/i B in V^Sgj/6Z-SiCNR) in the magnetic moments of these atoms; 
the third C atom surrounding (labeled 7 in Fig.[T]) has magnetic moment of about 0.80/xb. 
On the other hand, Fig. [8(b) shows that bond reconstructions in the vicinity of Vgj, which 
were observed in V^V^/GZ-SiCNR, do not occur when N is doped at the vacant C 8 site. Such 
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a process results in a increase in the local magnetic moments of the C atoms surrounding Vg; 
(we found magnetic moments of 0.17, 0.87 and 0.90/z.b in these atoms) and, consequently, 
in a increase in the magnetization of the system (Vg-Vc/GZ-SiCNR has magnetic moment of 
1.34/xe, while V^A^/GZ-SiCNR has magnetic moment of 3.06yU#). 

IV. CONCLUSIONS 

Ab-initio calculations were performed to investigate the magnetism of pristine and de- 
fected Z-SiCNRs. Single (Vsi and Vc) and double (VsiVsi and VsiVc) vacancies were con- 
sidered to probe the electronic structures of Z-SiCNRs. Our results indicate that these 
native defects can induce large magnetic moments in Z-SiCNRs. While the half-metallic 
character of the pristine Z-SiCNR is maintained in the presence of a single Vsi, we found 
that a single V c leads to a transition from half-metallic to metallic behavior in Z-SiCNRs, 
and double VsiVsi and V^Vc can transform half-metallic Z-SiCNRs into semiconducting sys- 
tems. Such electronic phase transitions were discussed in terms of hole doping. Single Vsi 
induced the largest magnetic moment among the defects investigated, but this defect has 
higher formation energy and is energetically unfavorable when compared with the single Vc 
and double VsiVc- The ferromagnetic ground state was shown to be more stable than the 
antiferromagnetic state; therefore, room-temperature ferromagnetism was also speculated 
in defected Z-SiCNRs. To reduce the defect formation energy of single V S i and realize it 
experimentally, the interactions of substitutional B and N impurities with this native defect 
were also studied. We found that a B atom substituting a Si atom, as well as a N atom 
substituting a C atom in the presence of vacancies, leads to a considerable reduction of the 
defect formation energy of Vsi- Therefore, we believe that light elements are beneficial for 
the realization of room-temperature magnetism in defective Z-SiCNRs. 
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